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Abstract 

Delays associated with elementary interaction processes are investigated. 
^. ■ The case of broad resonances is discussed in the context of reaction sim- 

O . ulations. 

°: 

q ■ 1 Introduction 

t^ ■ 

q\ . This study has been first motivated by problems within the transport theory 

for low-energy heavy-ion reactions. The common microscopic models for the 
reactions include the cascade and the Boltzmann-equation models. In the 
cascade model, a reaction is represented as a superposition of independent 
NN collisions. The Boltzmann-equations models account for the mean-field 
effects, in addition to collisions, with particles moving along curved rather 
than straight-line trajectories in-between the collisions. Within the models 
there is a flexibility in including different produced particles; the mean-field 
parameters may, principally, be constrained using flow data. 

In simulating the reactions, it is necessary to decide how to simulate the 
elementary scattering processes. Some possibilities are indicated in Fig. |l|. 
The most common prescription is the one where particles change their mo- 
menta at the distance of the closest approach. This constitutes the so-called 
impact scattering. Variations exist, correlating the spatial separation between 
the particles and the momentum transfer, as in the attractive (qr < 0) and re- 
pulsive (qr > 0) impact scattering. At low energies cross sections are isotropic 
and the choices include the billiard-ball scattering. Finally, the billiard-ball 
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Fig. 1: Different scattering prescriptions in transport models: impact-parameter, 
billiard-ball, and mixed scattering. 

scattering may be preferred at low impact parameters, and the impact scat- 
tering at high. This gives a so-called mixed scattering prescription. 

While it is difficult to assess which scattering prescription should be fa- 
vored for interactions, the results of reaction simulations exhibit surprising 
sensitivity to the prescriptions. Thus, in [ij variations by more than a factor 
of 2 in the rise of density above normal were found when using the billiard-ball 
and the impact-parameter scattering for nucleons in a central Ne + U reac- 
tion. In || the average flow angles were found to vary by nearly a factor of 2 
depending whether the impact or the repulsive impact scattering was used in 
a U + U collision simulation. Concerning high energies, Kahana et al. found 
huge flow effects in the simulations of Au + Au collisions at 10.7 GeV/nucleon, 
when utilizing the billiard-ball scattering. Herrmann and Bertsch [|| consid- 
ered an expanding pion gas in connection with the ultrarelativistic collisions 
and obtained large differences in source radii, such as studied through inten- 
sity interferometry, depending whether pions followed attractive or repulsive 
impact scattering. 

These and other experiences prompted an investigation |J to limit the 
short-term nature of elementary processes and establish possible links between 
the scattering prescriptions and the macroscopic properties of a system. 

2 Wave- Packet Interactions 

To assess what happens in scattering we consider an experiment where a target 
is located at the center of a scattering chamber of radius R, see Fig. ^. A wave- 
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Fig. 2: Scattering chamber of radius R. 



packet is directed, towards the target, asymptotically given by 



m = JdEg(E) L°(r,t) + f(0,E) 
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The wave-packet has a transverse size B and is constructed in such a way that 
at the time t = it is centered at the target. 

We examine the average time spent by the packet within the chamber 
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and compare it to the the time t® o1 = 2R/v spent within the chamber in the 
absence of interaction, r vo i = r® ol + At. In (|j), j is current density 



J 



— 0* d 
2fii 



(5) 



and N is integrated incident flux through the surface of the chamber. The prod- 
uct j ■ r in (U) is positive when a particle exits the volume and negative when 
a particle enters. 

For a large volume there is no interference between incoming and outgoing 
waves 



'vol 'out 'it 



(6) 



and the incoming time is found simply from (0) as Tj n = —R/v. The outgoing 
waves consist of scattered and unscattered portions of the wave packet and the 
outgoing current has contributions from the scattered wave, unscattered wave, 
and the interference between the two. 

To appreciate the importance of interference one may consider the differ- 
ence between the net number of particles going out from and coming into the 
chamber 

R 2 Jdt(Klj(R,n,t)-T = N out -N in . (7) 

By integrating over the back of the chamber, the current coming in (g(E) ~ 
6(E)) is found equal to 
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The current going out consists of scattered and forward parts, 

N out = N S + N f , (9) 

and 

N s = JdQ^ = JdQv\f\ 2 = va. (10) 

If the forward part is calculated from the unscattered packet only, the current 
coming out in the forward direction is found to be the same as that com- 
ing in (H). Thus, the current coming out through wide angles ( [TO]) appears 
to be born within the chamber. The solution lies in a careful treatment of 
the interference between the unscattered portion of the packet and the scat- 
tered wave. In (|5]), the interference gives rise to terms proportional to e lk ( z ~ r > . 
For finite angles the terms vanish upon an averaging over a small macroscopic 
range of R but that is not the case in the forward direction. On integrating 
around the forward direction with the inclusion of interference term, one finds 

N f = T2\ E( 2 ^ + 1 )- 4fcIm /(°) = v {^B 2 -a)=N m -N s . (11) 






Due to a = -^Im/(0), the interference terms precisely correct for the flux 
moved to finite angles. 

Returning to the times, at finite angles there is no interference and the 
average exit time may be computed as a function of scattering angle as 
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M ' dNjdfl v 2i\f\ 2 \ f dE f dEJ v dE ' 

(12) 



The information on the delay Ar s , relative to the case without interaction, is 
contained in the phase of the scattering amplitude. For scattering in only one 
partial wave the delay time is 

Ar. = ^. (13) 



Calculation of the average forward exit delay time yields 
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The delay time in the forward direction involves the real part of the forward 
scattering amplitude. 

3 Forward Delay and Mean Field 

One may ask whether any effect of the delays in elementary processes is ac- 
counted for in the reaction simulations. In the impulse approximation, the real 
part of optical nuclear potential may be expressed in terms of the real part of 
the forward scattering amplitude, 

TV 2ttv 
U = nReT(0) = - v -^Ref(0), (15) 

indicating a possible connection to the forward delay time. Note that, when 
density varies in space, so does the optical potential. 

To further explore the possible connection between the optical potential and 
the forward delay, let us consider a situation where we construct an optical 
potential landscape associated with some distribution of scatterers in space. 
One way to do it is to surround each scatterer with a well of radius R and depth 
Vo = ReT(0)/V, where V is the volume of a well, and to increase R until the 
wells merge, see Fig. [| Let a particle be incident on a set of the scatterers at 
a low density. The average time spent by the incident particle within the region 
of scatterers will be affected by the optical potential. The time of traversal 
through a single well is r = s/v, where s is the length of trajectory within 
the well and v is group velocity. When the optical potential is put to zero, 
the average length of a trajectory becomes just the ratio of the volume of the 
well to the transverse area, s = V/A = | tc R 3 /tc R 2 = | R. The average time 




Fig. 3: Particle incident on potential wells. 

within the well is r = s/v. When the potential is switched on, the time within 
the well changes compared to the case without potential, due to a change in 
the length of trajectory caused by refraction, and due to the change in the 
group velocity, 
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The above manipulation shows that the change in the average time for a single 
well is identical to the forward delay for a packet of transverse size equal to 
the size of the well, X^(2£ + 1) « k 2 R 2 . Thus, the forward delay time and the 
mean field are, indeed, directly related. 

4 Delays and Pressure 

Pressure for an enclosed system may be defined as the momentum delivered to 
the walls per unit area and per unit time. Given a dependence of the pressure 
on density and temperature one can determine other thermodynamic quantities 
for a system. In the mean-field approximation, at low density, the correction 
to the ideal-gas pressure from interactions can be, conversely, obtained from 
the energy density. Simple manipulations allow to express the correction to 
the pressure in terms of the forward delay time, 
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Given that particles stay shorter or longer within the interaction region, com- 
pared to the free-flight time, they knock more or less often against the en- 
closure walls and the pressure increases or decreases, respectively, due to the 
interactions. Depending on the sign of the delay time the correction due to the 
interactions may be further interpreted in terms of excluded volume (At < 0) 
or a reduction in the number of degrees of freedom (At > 0). 

Yet the delays associated with scattering should also affect the frequency of 
knocking against the walls. To the second-order in density the pressure is given 
completely by the Beth-Uhlenbeck formula. Manipulations of the second-order 
term in the formula separate out the contributions to the pressure from the 
forward delay time or mean field and from scattering, 
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The scattering can affect thermodynamic properties of a system just as the 
mean field can. Yet, the nuclear equation of state in connection with such 
reaction simulations as discussed in the introduction is usually considered in 
terms of the mean field alone. It is apparent from ( |T8"D that, if one aims at 
describing solely properly the thermodynamic properties, then, to the second 
order in density, one can entirely absorb the effects of scattering into the 
forward delay or the mean field or, alternatively, the effects of the mean field 
into scattering. The first is accomplished by using, in place of the forward 
delay Atj, the time 

» , . kR 2 r ,_ da . , . 

At =At ; + dn — AT s , (19) 



or by making a replacement in the mean field 
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The effects of the mean field may be absorbed into scattering by using, in place 
of AT,, the time 

Ar * , = ATs + ^ Ar/ - (21) 

At a general level, the relation of thermodynamic properties to the delays 
may be understood in terms of ergodicity. The system spends time in a certain 
phase-space region, that is proportional to the density of states in that region. 
Changed time corresponds to the changed density of states and thus to changed 
thermodynamic properties. 

4.1 Resonance Interactions 

One of the important examples of interactions to consider is that of a Breit- 
Wigner resonance. The phase shift for the resonance is given by 
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In transport simulations the resonances are given, intuitively, a lifetime At s = 
1/r. From ( |22"D the actual scattering delay is 

ds r 

s = dE = 2[{E-E R y + ry4] ' ( } 

At E = Er, the scattering time is twice as large as the naively expected and 
it tends to zero as \E — Er\ gets large, see Fig. (|. The forward delay time, 
normalized in such a way as if it were to be put into scattering, is 

nB 2 a , A (E - E R ) 2 - r 2 /4 , . 

At/X _ = At ;- At , = _|_^__^_ . (24 ) 

This time is negative in the vicinity of resonance. That corresponds, in par- 
ticular, to the known increase in the group velocity for light near resonances 
in dielectrics. 



The sum of the times fl23| ) and ([24] ) is Ar' s = 1/r. Thus the naive scattering 
time accounts for both the scattering and the mean field. When employing such 
a time in a simulation it is inappropriate to include separately the effect of the 
specific resonance onto mean field. For narrow resonances, under averaging 
over energy, the time At, averages itself out to 1/r, as the mean field averages 
out to zero. However, the averaging would not be justified for broad resonances. 

One of the more important broad resonances in high-energy nuclear physics 
is the 7riV A-resonance. The width for this resonance is comparable to the 




Fig. 4: Time delays in the case of a Breit-Wigner resonance as a function of 
energy from the resonance divided by the width. The solid line represents the time 
for the scattered wave At,. The short-dashed line represents the time At' s which 
is ergodically consistent when the delay for forward wave or the mean field are 
neglected. Finally, the long-dashed line shows the time delay for the forward wave 
divided by the fraction of the incoming wave that is scattered, At' s — At s . 



energy above the threshold and, as such, the width exhibits a strong energy 
dependence. In the immediate vicinity of the threshold, the width behaves as 
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The use of the time 1/T(E) in simulations leads to very long-lived A's close 
to the threshold, which must be unphysical. Intuitively, the A's far from the 
center of the resonance are expected to be virtual. Figure [| shows the delay 
times At s and At' s — At s computed from the measured irN phase shifts. It is 
apparent that the scattering delay time Ar s has a quite regular behavior with 
energy. It is At' s — Ar s , within At' s ~ 1/r, that strongly diverges close to the 
threshold for the A resonance production. It then follows that the physically 
unacceptable divergence of At' s in simulations is associated with forcing the 
mean field effects onto scattering when actually the scattering cross section 
declines, a n N->A oc (E — m n — Mat) 1 / 2 . 

With the above, the solution to the dilemma with diverging times appears 
straightforward. The scattering and forward delays should simply be treated 
separately. Yet the delays in a simulation accomplish various goals. One goal is 
the production of the proper overall density of states and of the corresponding 
proper thermodynamic properties in equilibrium. Another goal is the produc- 
tion of the proper number of resonances which, when decaying, may contribute 
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Fig. 5: Time delays for a n + p system as a function of cm. kinetic energy, computed 
using measured phase shifts. The solid line represents the time delay for the scattered 
wave averaged over angles and spin directions, Ar s . The dashed line represents the 
forward time delay averaged over spin directions and divided by the fraction of the 
incoming wave that is, on the average, scattered, A.t' s — At,. 



to the energetic photons or lepton pairs, etc. In equilibrium, the number of A 
resonances is given in terms of the resonance spectral function and the use of 
the time At' s = 1/T in scattering, during which a irN pair converts into a A, 
meets the requirement of the proper number of resonances, see below. One 
may expect that, when replacing At' s by Ar s and Arj, and when converting 
ttN pairs into the A's during these times, the requirement will still be met, 
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However, the problem which arises, is that the forward delay time for a res- 
onance (|2~3| ) is not positive definite, and one cannot simulate a negative time 
for the conversion of a nN pair into a A. A modification, relying on the pos- 
sibility of moving the mean-field effects into scattering and back, can render 
the scattering and forward times both positive and free from any singularity, 



e.g. with 
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For the forward processes the A should decay back into the original pion and 
nucleon with their original momenta. The times ([27|) are consistent with the 
ergodic constraint, up to terms proportional to dT /dm and dm A /dm, and 
with the proper number of A's, as related to the number of tcN pairs. 

5 Schemes in Simulations 

The sensitivity of the reaction simulation to prescriptions in scattering, dis- 
cussed in the introduction, can be understood in terms of scattering delays 
and changes made to the nuclear equation of state. The billiard-ball scat- 
tering gives an average delay time equal to the negative of the time within 
hard-core in a free passage, 

1 |vrci 3 Ad 

Ar. = — S-— = ---. (28) 
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The insertion of this time delay into the pressure ( jig) yields the repulsive 
excluded volume correction to the pressure, valid to within the second order 
in density, 

P = nT--n 2 T(av Ar s ) =nT + n 2 T — . (29) 

The sign of the time delay and of the correction to pressure may be inverted, 
making the correction to pressure attractive, by replacing the hard core in 
the scattering by a spherical shell of radius d open in the direction of motion, 
(essentially replacing the convex by the concave mirror). 

The repulsive impact-parameter scattering yields the average time delay 
and the correction to pressure nearly as large as the billiard-ball scattering, 

Ar. = -£-, (30) 

3 v 

i.e. 7r/4 of (p8|). The attractive impact-scattering changes the sign in ( |30| ) and 
in the correction to pressure. Finally, for the impact scattering the average 
time-delay is zero and the correction to pressure from scattering vanishes, 
within the second order in density. Given required delay times obtained from 
phase shifts 0, such that |Ar a | < (2/3) d (sin6>) = (Ar s ) max , these times may 



be generated making a fraction v of all scatterings repulsive and a fraction 
1 — v attractive, 
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6 Conclusions 



The two- or more-body processes hardly ever take place instantaneously. In an- 
alyzing the wave-packet scattering we found that, in the consequence of inter- 
action, both the forward and the scattered waves get delayed. The forward 
delays are accounted for in the mean field in reaction simulations. Both types 
of delay, though, through ergodicity, affect the density of states in a system 
and the system thermodynamic properties. To the extent that only the density 
of states matters, then, to the lowest order in spatial density, the same density 
of states and the same thermodynamic properties may be obtained by forc- 
ing all delays onto scattering or onto the forward direction. Some difficulties, 
though, may arise in specific situations. Thus, the unphysically long lifetimes 
of resonances close to threshold in simulations are specifically associated with 
the attempt to force the mean-field effects onto scattering. 
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